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We investigate decoherence and equilibration in the experimentally relevant situation of weak 
coupling to an environment. We consider small subsystems of large, closed quantum systems that 
evolve according to the von Neumann equation. Without approximations and without making any 
special assumptions on the form of the interaction we prove that, for almost all initial states and 
almost all times, the off-diagonal elements of the density matrix of the subsystem in the eigenbasis of 
its local Hamiltonian must be small, whenever the energy difference of the corresponding eigenstates 
is larger than the interaction energy. This proves that decoherence with respect to the local energy 
eigenbasis is a natural property of weakly interacting quantum systems. 

PACS numbers: 05.30.-d, 03.65.-w, 03.65.Yz 

Keywords: open quantum systems, equilibration, non-Markovian dynamics, einselection, pointer states 



INTRODUCTION 

Quantum Mechanics claims to be a fundamental the- 
ory. As such, it should be able to provide us with a 
microscopic explanation for all phenomena we observe 
in macroscopic systems, including irreversible processes 
such as thermalization. But its unitary time development 
seems to be incompatible with irreversibility, leading to 
an apparent contradiction between Quantum Mechanics 
and Thermodynamics jl8j . 

To explain irreversible processes and to overcome the 
problem of the emergence of classicality many authors 
have suggested to modify Quantum Theory. For exam- 
ple, by adding non-linear terms to the von Neumann 
equation, or by postulating a periodic spontaneous col- 
lapse of the wave function [T] ^2] . Others have considered 
Markovian, nonunitary, time evolution [2,i3| which can be 
thought of as resulting from an interaction with a memo- 
ryless bath and it has been shown that system bath mod- 
els that evolve under a special type of Hamiltonian tend 
to evolve into states that are classical superpositions of so 
called pointer states [4| — a phenomenon called einselec- 
tion. These approaches, which are subsumed under the 
term decoherence theory, are able to reproduce many of 
the features of dissipative systems and are undoubtedly 
very valuable for applications. 

But, in face of the enormous success of standard Quan- 
tum Mechanics in explaining microscopic phenomena and 
the existence of macroscopic quantum systems on the 
one hand and the broad applicability of Statistical Me- 
chanics and Thermodynamics on the other, we feel that 
neither a modification of Quantum Theory, nor consid- 
erations restricted to special Hamiltonians can provide a 
satisfactory explanation for the classical, statistical, and 
thermodynamic behavior of the macroscopic world. Re- 
cently there has been remarkable progress in explaining 
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macroscopic, seemingly irreversible behavior from stan- 
dard Quantum Mechanics. It has been shown that it is 
possible to explain the phenomenon of equilibration and 
irreversibility [5i-(7j and to justify the applicability of the 
canonical and microcanonical ensemble |5] without added 
randomness (i.e., without assuming the existence of al- 
ready equilibrated and thermalized baths) and ensemble 
averages, from nothing but pure Quantum Mechanics and 
the randomness due to entanglement with the environ- 
ment (see also [9j and the references therein). 

We make use of the results obtained in these papers 
and connect this approach with the research on decoher- 
ence. We consider the case of decoherence due to weak 
interaction with an environment. A weak coupling to 
an environment exists in practically all situations. This 
case is thus of particular interest for developing a better 
understating of the foundations of Statistical Mechan- 
ics and Thermodynamics, for applications in quantum 
information processing and quantum computing and for 
experiments on environment- assisted entanglement cre- 
ation (see [10] and the references therein). Our main 
result is that decoherence with respect to the local en- 
ergy eigenbasis is a natural property of weakly coupled 
systems. 



SETUP AND NOTATION 

We consider arbitrary quantum systems that can be 
described using a Hilbert space H of finite dimension d 
and that can be divided into two parts, which we will call 
the bath B and the subsystem S. For infinite dimensional 
systems it is often possible to find an effective description 
in a finite dimensional Hilbert space by introducing a 
high energy cut-off. Moreover, it was demonstrated in 
[TTj that many of the phenomena that can be rigorously 
proven in the finite dimensional case also occur in infinite 
dimensional systems. We thus believe that the restriction 
to finite dimensions as mainly a technicality. 

We use the terms bath and subsystem because in the 
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end we will be interested in situations where the dimen- 
sion ds of the Hilbert space of the bath "Hs is much 
larger than the dimension ds of the Hilbert space Hs of 
the subsystem, besides that S and B are two completely 
arbitrary quantum systems. 

We assume that the Hamiltonian of the joint sys- 
tem has non- degenerate energy gaps. This assumption al- 
ready appears in the work of von Neumann |12j and later 
in [51 in] and means that for any four energy eigenvalues 
Ek,Ei,Em, En equality of the gaps E^ - Ei ^ E„^ - E^ 
implies that either k = I and m = n ot k = m and 
I — n. It shall be emphasized that this is an extremely 
weak assumption and with some additional effort, it can 
be replaced by an even weaker one that allows degen- 
eracies of the energy levels [7]. Every Hamiltonian be- 
comes non-degenerate by adding an arbitrary small ran- 
dom perturbation; therefore the Hamiltonians of macro- 
scopic systems can be expected to satisfy this constraint. 
The physical implication of the above assumption is that 
the Hamiltonian is fully interactive in the sense that there 
exists no partition of the system into two non interacting 
subsystems. For fully interactive systems our results are 
robust against the existence of some degeneracies in the 
energy gaps. How non-degenerate the energy spectrum 
is influences the equilibration and decoherence times. 

We use p for density matrices of possibly mixed states 
and tp if the state is pure. All states are assumed to be 
normalized Tr[p] = 1. Their reduced states on the bath 
and subsystem are denoted using superscript letters like 
in — Ti' s[p] and p^ = TrB[p]- We write the trace 
norm of a density matrix p as ||jo||i = Tr[^ p^ p] — Tr 
and the trace distance as 

V{p,<j)^^\\p-a\U. (1) 

We denote the operator norm of a Hermitian operator A 
acting on some Hilbert space H by 

\\A\\oo^ max Tr[^V], (2) 

where 7^i('H) is the set of rank one projectors on H. We 
use the letter w to denote the time average of time de- 
pendent states pt 

^ = {Pt)t = lim - / ptdt. (3) 

r-)-oo r Jq 

EQUILIBRATION 

In a time reversal invariant theory equilibration in the 
usual sense is impossible. We therefore use an extended 
notion of equilibration and say that a system is in equi- 
librium when its density matrix stays close to some state 
for almost all times, and say that it evolves toward equi- 
librium if it approaches such a state and then stays close 
to it if started in a state far from equilibriimi. 



Recently, it has been shown that under the above as- 
sumptions and whenever the initial state has a high ef- 
fective dimension d^^{uj) = 1/Tr[a;^], uj = {pt)t every 
small subsystem with ds ^ equilibrates in this 

extended sense: 

Theorem 1 ([6 ). Consider any pure state -04 evolving 
under a Hamiltonian with non- degenerate energy gaps. 
Then the average distance between pf = Tr^ il>t and its 
time average ui^ = {pf)t is hounded by 

The second important result of [6] is that the effec- 
tive dimension d'^^{uj), which is a measure for how many 
energy eigenstates contribute significantly to the initial 
state, is large for almost all pure states drawn according 
to the unitary in invariant Haar measure: 

Theorem 2 ([6 ). (i) The average effective dimension 
{d'^^ {oj)) , where the average is computed over uniformly 
random pure initial states ■00 G Vi{'Hr) chosen from a 
subspace Hu of dimension dn, is such that 

(d^«(c.)>^„ > ^. (5) 

(ii) For a random pure initial state t/jq G Vi{'Hr), the 
probability that d°^{uj) is smaller than d^jA is exponen- 
tially small, namely, 

Pr|d<^«(u;)<^|<2e-^^ (6) 

with a constant C = \n{2)'^ / {72 tt^) . 

The Haar measure used in the above theorem is some- 
times criticized for being unphysical. As the results pre- 
sented herein depend crucially on d'^^ being large it is 
therefore worth saying a few words about why we believe 
in their physical significance despite the criticism con- 
cerning the Haar measure: first theorem [5] is a measure 
theoretic result, it is not to be misunderstood as state- 
ment about states drawn from an actual physical ensem- 
ble. The bound on the probability to get a state with 
a low effective dimension drops off exponentially. This 
raises the hope that the result does not depend on the 
details of the measure from which the states are sampled 
and that similar statements can be proven for other non- 
singular measures. Second, theorem [2] is a very strong 
statement and what is actually needed in the following 
is much weaker, namely that d'^^{u!) is much larger than 
some low, fixed power of ds, which can be as low as 4 
or 8 for a single Qubit. It seems to be unreasonable to 
assume that the quantum state of a macroscopic object 
is composed of only that few energy eigenstates. 
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SPEED OF FLUCTUATIONS AROUND 
EQUILIBRIUM 

Knowing that, under suitable conditions, subsystems 
of large quantum mechanical systems will equilibrate, it 
is natural to ask: how fast will the fluctuations around 
the equilibrium state typically be? This question was 
investigated very recently in [7] . 

The first step is to introduce a meaningful notion of 
speed. This is achieved by defining the time derivative [7j 



vsit) 



lim 

St^O 



St 2 



dpi 



dt 



with 



dpi 
dt 



i TrB[p4, Jf] 



(7) 



(8) 



As the choice of the origin of the energy scale does not 
influence the speed, it is convenient to split up the Hamil- 
tonian of the system in a part J^q proportional to the 
identity and the traceless operators J^s^ b and sb 
as follows: 

= Jfo + ^S'»l + l«)^i3+^Si3 (9) 

Using a result from [5] it is shown in [7] that: 

Theorem 3 ([?])• For every initial state po of a com- 
posite system evolving under a Hamiltonian of the form 
([9]) and with non- degenerate energy gaps, it holds that: 



{vsit))t < |l^s®l + ^ 



SB II00I 



dl 



d<=ff(u;) 



(10) 



As we have argued above, we are convinced that the 
effective dimension d'^^{uj) is typically very large in re- 
alistic thermodynamic systems. In particular, as all di- 
mensions grow exponentially with the number of con- 
stituents of the system it will usually be much larger 
than any fixed power of ds- Therefore, the speed of the 
subsystem will, most of the time, be much smaller than 
II J^s €5l -I- J^SB II 00, which in turn can be expected to 
grow at most polynomial with the number of constituents 
of the subsystem and is the natural unit in which the 
speed of is to be measured [7] . 



states \p), invariant [H]. If this is the case, the Hamilto- 
nian and the time evolution operator have the form 



p 



(p) 
t 1 



(11) 

(12) 



where Ui 



ip) 



are arbitrary Her- 
mitian matrices. One finds that the subsystem state of 
an initial product state of the form po = Pq ipQ , where 
the state of the bath can be assumed to be pure without 
loss of generality, evolves into 



Pt 



(13) 



Under the evolution induced by such a Hamiltonian the 
diagonal entries of Pq, when expressed in the pointer ba- 
sis, remain unchanged while the off-diagonal entries are 



suppressed by a factor of {iJjq \ u'f ' 



r(p) 



t 



< 1. The 



actual time development of the {iPqIU^^ U^^^i^q) de- 
pends on the explicit model under consideration, but for 
many models they have been found to decrease rapidly 
over short time scales [HI [H [H]. If some of the Jf^^^ 
lead to an identical time development for the chosen ini- 
tial bath state there exist subspaces of T-Ls in which co- 
herence is preserved. 

Note that, the diagonal entries, which survive the de- 
coherence, are completely determined by p^ and do not 
depend on the initial state of the bath ip^ at all. The di- 
rect opposite situation is the thermodynamic case where 
the final state is completely determined by the properties 
of the bath. Most realistic situations surely lie between 
these two extremes. 



EQUILIBRATION AND EINSELECTION 

Using the results discussed in the preceding sections it 
is possible to get rid of the quite limiting assumptions on 
the Hamiltonian and to show that einselection, i.e. de- 
coherence with respect to a fixed basis, naturally occurs 
in situations where the interaction Hamiltonian M' sb is 
weak. 

According to [7) the velocity can be written as 



EINSELECTION IN A NUTSHELL 

The term einselection^ which stands for environment- 
induced super selection, is due to Zurek [H [13] . Einselec- 
tion is known to occur in situations where the Hamilto- 
nian of the composite system leaves a certain orthonor- 
mal basis of the subsystem, spanned by so called pointer 



J S s 

k=l 



(14) 



where the dg operators Ck form a Hermitian orthonormal 
basis for Hs such that Tr[efc e/] = Ski and 

Ck{t) = Tr [p{t) i [^5 ®1 + y^SB,ek ® 1]] . (15) 
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The velocity depends on b only implicitly through the 
trajectory pt, but for an arbitrary fixed state p the ve- 
locity is solely determined by ^5 and J^sb- 



Proof. Using the inverse triangle inequality and ( 16 1 we 
see that 



lit 



i[p'\j^s]\\i-\\iTTB[p,J^SB]\\i\ < 



-i[/,Jf's]+iTrs[p,^SB] 



(16) 



dp'= 



dt 



(22) 



Now if J^sB is much weaker than (16) is domi- 

nated by the first term. Consequently, the system can 
only become slow when [p^ , J^s] is small. To see when 
this happens we first establish a general lower bound on 
the norm of commutators between states and arbitrary 
Hermitian matrices: 

Lemma 1. Let p be a normalized state and A a Her- 
mitian observable with eigenvalues Uk and eigenvectors 
lofc), then 



For to become small the norms of the two 

commutators must be approximately equal. Applying 
lemma [1] to the norm of the first commutator yields: 

II i[/,^s]||i > 2 max|£;f ~ i?f | |pfj (23) 

The norm of the second commutator can be upper 
bounded, using the well-known fact that the trace norm 
of traceless, Hermitian matrices is non-increasing un- 
der completely positive, Hermitian, trace-non-increasing 
maps as follows: 



',A]||i = ||i[p,A]||i >2 max ^|afc-a,||pfei| (17) \\i Ti-b[p,J'^sb]\\i < \\[p, sb]\\i < 2\\ J^^sb \ 



> 2max|afc - a;| \pki\. 

kl 



(18) 



(24) 



□ 



where the maximization is performed over all decompo- 
sitions of the index set {1, . . . , ds} into non overlap- 
ping pairs (k, I) over which the sum is performed and 
Pki = {ak\p\ai)- 

Proof. The equality is trivial. For all traceless, 
Hermitian, bounded operators B on some finite di- 
mensional Hilbert space T-L it holds that ||i?||i — 
2 maxng-p(-H) Tr[Hi?], where P(H) is the set of all pro- 
jectors on T-L and the maximum is obtained when H is the 
projector onto the positive subspace of B. By expand- 
ing p in the eigenbasis of A, using the above equality for 
B = [p, A\ and considering all sums of orthogonal rank 
one projectors Ti^i of the form 



H 



hi 



kfc;)(7rfci| Itth) 



1 



72 



(lafc)-fe'^-la,)), (19) 



where 4>ki are phase factors, one easily verifies (17). The 
second inequality is trivial. □ 

We can now prove the main result of this paper: 

Theorem 4. Consider a physical system evolving under 
a Hamiltonian of the form given in ^ and with non- 
degenerate energy gaps. All reduced states p^ satisfy 



SB 



dp^ 



dt 



> max V|i?,^ -i?f||p^;| 



{(k.i)} 



> max \Ei 
~ kl ' 



(20) 

EfWpil (21) 



where pf^ = {E^\p^\Ef) and and \E^) are the eigen- 
values and eigenstates of J^s- 



The assertion of theoremj4] is almost intuitively clear, 
but combined with theorem[3lit allows to draw the follow- 
ing important conclusion: Whenever d°^{uj) is large the 
subsystem is slow most of the time and if this is the case 
coherent superpositions of eigenstates of J^s with eigen- 
value differences that are much larger than || J^sb \\oo 
can not contribute significantly to the state of the sub- 
system. That is, the corresponding off-diagonal elements 
of the reduced state pf in the Jifs eigenbasis must be 
small. A similar behavior was observed for a specific 
model in [TS] |T5] . 

Without using any approximations we have shown 
that coherence can only be retained between eigenstates 
of J^fs whose energy difference is small compared to 
II '^SB lloo- This statement remains meaningful even 
when the subsystem is large and its energy spectrum 
thus very dense. Theorem |4] then still implies that coher- 
ent superpositions of eigenstates with far apart energies 
(sometimes called Schodinger cat states) must decohere. 
If the subsystem is small and the interaction Hamilto- 
nian weak compared to the energy gaps of the subsystem 
Hamiltonian it implies an even stronger statement. The 
state of the subsystem must then, most of the time, be 
approximately diagonal in the eigenbasis of ^s- 



CONCLUSIONS 

We have shown that quantum systems which interact 
weakly with the environment tend to evolve into con- 
vex combinations of energy eigenstates of their Hamil- 
tonian. This result is obtained without making any ap- 
proximations and without neglecting memory effects in 
the bath. No assumptions on the details of the inter- 
action are made, as opposed to the classical einselection 



5 



mechanism due to Zurek. In particular, we do not need to 
assume that the interaction with the environment leaves 
a certain set of pure pointer states invariant. This proves 
that decoherence with respect to a fixed basis is a very 
natural property of weakly interacting quantum systems. 
Due to the generality of our approach we cannot say much 
about the time scales on which decoherence happens. For 
this, specific models must be considered [TH [17] . Our re- 
sult establishes a link between decoherence theory and 
the recent research on equilibration and the foundations 
of Statistical Mechanics [HI [7] . 

Decoherence in the energy eigenbasis is observed in 
many situations where the local Hamiltonian is much 
stronger than the interaction. A well-known example 
is electronic excitations of gases at moderate tempera- 
ture. The energy gaps between the ground state and the 
first few excited states are typically much larger than 
the thermal energy. The dynamics of such systems is 
successfully described using transition rates between en- 
ergy eigenstates. Ultimately theorem |4] explains why this 
is eligible. 
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